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We examine spectral equilibration of quantum chaotic spectra to universal statistics, in the context 
of the Brownian motion model. Two competing time scales, proportional and inversely proportional 
to the classical relaxation time, jointly govern the equilibration process. Multiplicity of quantum 
systems having the same semiclassical limit is not sufficient to obtain equilibration of any spec- 
tral modes in two-dimensional systems, while in three-dimensional systems equilibration for some 
spectral modes is possible if the classical relaxation rate is slow. Connections are made with upper 
bounds on semiclassical accuracy and with fidelity decay in the presence of a weak perturbation. 
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I. INTRODUCTION 



CN , Random Matrix Theory (RMT) was introduced by Wigner in the 1950s as a model for describing the universal 
' spectral behavior of complex many-body systems such as compound nuclei. It was conjectured that a typical spectrum 
of a complex system displays statistical properties similar to those of a Hamiltonian chosen at random from a basis- 
independent ensemble of random matrices in the same symmetry class. Subsequently, Dyson showed that the statistical 
predictions of RMT could be reproduced within a Brownian motion model, where Brownian motion of the energy 
levels results from a stochastic perturbation of some initial, possibly nonrandom, Hamiltonian Bohigas, Giannoni, 
^ ' and Schmit argued that the same universal behavior of the spectrum occurs generically even in single-particle systems, 
1^ . as long as the underlying classical dynamics is chaotic 0. 

^ ^ ' Wilkinson used the Brownian motion model to explain the connection between quantum chaotic systems and RMT 
spectra, using the "quantization ambiguity" (i.e, multiplicity of quantum systems having the same semiclassical limit) 
' as the source of perturbation causing Brownian motion of the energy levels and eventual equilibration to RMT spectral 
J> \ statistics . More recently, the Brownian motion approach has been extended to study the evolution of eigenstates 
^^O ■ and to relate the smooth and oscillatory parts of the spectral correlation function in diffusive and ballistic chaotic 
' systems 0- 

. It well known that not all spectral modes of a typical quantum chaotic spectrum obey universal statistics. It is 
therefore of interest in the context of the Brownian motion model to investigate which modes are equilibrated in a 
, given dynamical system, for a given perturbation strength (such as that arising from the quantization ambiguity), as 
a function of degree of chaos in the system, and as a function of effective h. Alternatively, one may ask about the 
j size of perturbation necessary to achieve equilibration to universal statistics for a given spectral mode; this question 
. ^ ■ is related to work by Zirnbauer on energy level corr elations in a quantum ensembre of ckssically identical maps % 
(— I \ as well as to the method developed by Gornyi and Mirlin for studying wave function correlations in ballistic systems 
by adding weak smooth disorder 

In the present work, we show that the Brownian motion model can be used successfully to predict the number of 
equilibrated spectral modes as a function of perturbation strength, chaoticity of the classical system dynamics, system 
%J^ • dimension, and effective h. The equilibration process is governed by two competing time scales, allowing the number of 
5^ , modes equilibrated to RMT either to increase or to decrease as the original system becomes more chaotic; equilibration 
is reduced for systems with very short or very long classical relaxation times and is maximized for systems with 
intermediate Lyapunov exponent. We show that in two-dimensional quantum systems no equilibration to universal 
behavior for any spectral modes can arise from the quantization ambiguity in the ?i ^ limit; in three dimensions 
some modes may be equilibrated depending on the classical relaxation time of the underlying classical dynamics; in 
four dimensions and higher, equilibration to universality will occur for at least some modes, but maximal equilibration 
requires the classical relaxation time to be much longer than the typical ballistic time scale of the dynamics. 

The paper is organized as follows. In Sec.m we review semiclassical estimates for matrix element variance, which are 
required as input for the Brownian motion model. In Sec. lIIII we discuss several example systems used for comparison 
between theory and numerics, and in Sec. II VI we review the Brownian motion model itself. Two competing conditions 
for equilibration of spectral modes are obtained in Sec. [V] and applied to equilibration of the nearest level spacing 
distribution in Sec. lVIl The implications for the effect of quantization ambiguity on spectral universality are developed 
in Sec. IVIll Finally, in Sec. IVIlH we demonstrate the relationship between results obtained in the Brownian motion 
model and recent findings on semiclassical accuracy Q and decay of quantum fidelity . 
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II. SEMICLASSICAL CALCULATION OF MATRIX ELEMENT VARIANCE 

This discussion follows 0,0, E|- Consider the matrix elements Bnn — (^n|-S|^'„) of an operator B, given the 
eigenstates |^„) of Hamiltonian H. It is assumed that the operator B is "independent" of H and has a well-defined 
classical limit B{q,p). We also assume that the phase-space average of the corresponding classical observable vanishes, 
so that the matrix elements fluctuate around zero, _B„„ = 0. The variance of the matrix elements at energy E is 



aliE) - S,L = HEr'J2 BIME - E^) , (1) 



where 5^{E) — i^TTrf) ^/^e ^ Z^** , v{E) — (5,,(-E — En) is the density of states, and • • • denotes averaging over a 
"suitable" ensemble. One may approximate 



BnnSe{E — En) 



(2) 



where e = V^rj is of the order of but smaller than the mean level spacing A — \/v{E). In the semiclassical limit, the 
expression in square brackets is approximated by a trace formula 

ZTTn ^ — ' 

p,r 

Here the sum is over all primitive periodic orbits p at energy E and their repetitions r, Tp is the primitive orbit 
period. Bp is the classical average of the observable B over the orbit, Dp^r = |det(Mp^j. — 1)|~^/^ is the square root 
of a classical focusing factor, Mp^r is the monodromy matrix of the orbit, Sp,r is the action, and /Zp.r is the Maslov 
index. Using (Im/)^ — iRe/*/ and neglecting repetitions, one obtains the following estimate for the variance, 

<^%{E) - ^^{2nh)''Y.T^\Bp\'Dle-^'^^/^' , (4) 

where as usual /3 = 1, 2 characterizes the presence or absence of time reversal symmetry, respectively. Now, assuming 
Bp to be real, and averaging over many periodic orbits of period close to Tp, 



Bl = r dti r dt2 B{ti)B{t2) r dt {BiO)B{t)) , (5) 

"'0 "'0 J-Tp/2 



where {B{0)B{t)) is the classical average of B{q,p)B{q{t),p{t)) over the energy hypersurface at energy E. Making 
use of the classical sum rule 

^ Tp ,f{Tp) Dl e-^^n/^' « 2 / dT f{T) e^^^^/^^ , (6) 
p Ja 

one finally obtains 

" 1^(2.^)"' r dT f'\t{Bit)Bme-^'^'/^^' 
P HE) Ja J-T/2 

2 1 



(3 TTu{E)n Jo 

with /e(<) = 1 — erf(2ei/?i). This may be approximated as 



dt{B{t)BiO))Mt), (7) 



where 



O /r2\ poo 

''BiE)^-^^jj^ / dtP{t)Mt), (8) 



, , (B(t)B(0)) 
Pit) = ^ (9) 
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is a classical correlation function. For a chaotic system, P{t) typically decays as a sum of exponentials [Tl|. The 
long-time behavior is then governed by the smallest exponent A, P{t) « 6e~^', where 6 is a classical constant of order 
unity. If the long-time behavior dominates the integral, one obtains in the limit e Xh 



2 b(B^ 



46 (B^) 



(3nv{E)hX (3XTh 



46 {B^) Tde 



(10) 



where Th — 27r?i/A = 2TThv{E) is the Heisenberg time, at which individual eigenstates are resolved, and Tdccay A 



is a classical correlation decay time. Since v (E) cx h^", the variance decays as 

al{E) cx fi'^-^X-^ 



(11) 



for small fi. This is consistent with Shnirelman's conjecture concerning the ergodicity of chaotic eigenstates. The result 
is independent of the smoothing function Se{E). In the limit of small A, the decay time Tdccay — A~^ becomes longer 
than the Heisenberg time Th = 2TThv{E), and A~^ in Eq. l(TT|l should be replaced with 2TTfiv{E). Then aj^iE) cx 
as expected for a quantum regular system. We can interpolate between the extreme chaotic and effectively regular 
regimes by taking Tdccay intermediate the period of the shortest orbit and the Heisenberg time. 



III. CLASSICALLY CHAOTIC MAPS 



As an example, consider the class of maps 



Pk+l = Pk - V'{qk) 
Qk+l ^ qk+ T'{pk+i] 



(12) 



on the torus [— tt, tt) x [— tt, tt) . Any such map may be thought of as arising from the periodically kicked Hamiltonian [l^ 

El 



The choice T{p) = vp'^/2 + K2vsmp and V{q) 
Here 



H{p,q,t)^-^T{p) + V{q) V ,5(t-jTkick). (13) 

J kick . 

j = -oo 

: —vq'^/2 — Kiv sing corresponds to a perturbed sawtooth map. 
Ki V cos qk 



Pk+i = Pk + vqk 
qk+i = (1 + v'^)qk +vpk + Kiv'^ cos qk + K2V cos pk+i 

and linearized motion is governed by the monodromy matrix 



M : 



1 



v{l- Ki sinqk) 



v{l- K2 sinpk+i) 1 + -y^ (1 - i^i sin(7fc)(l - K2 sinpfc 



(14) 



(15) 



In this case, strict hyperbolicity is guaranteed for v ^ 0, as long as \K\\ < 1 and \K2\ < 1. The parameters Ki and K2 
introduce nonlinearity and symmetry breaking into the dynamics. The classical Lyapunov exponent and correlation 
decay time Tdccay may easily be controlled by adjusting the parameter v. 

In Fig. ^ the classical autocorrelation function P{t) for the perturbed sawtooth map is shown for several values 
of the V parameter. Clearly, the long-time decay exponent A decreases with decreasing u, and empirically we find 
Xv ~ v^'^ for w <C 1. Note that the correlation decay rate A behaves very differently from the Lyapunov exponent, 
which scales linearly with v. 

When a discrete-time map is quantized on an iV— dimensional Hilbert space, we have h = 2t:/N, I'iE) = 
(iV/27r)2 Tkick, and Th = A^Tkick- Then Eq. ® becomes 



2 2(B^ 



and Eq. (|10|) takes the form 



13 N 



^ 00 



k=l 



2 2 6(^2 



/3AA^Tkick 



(16) 



(17) 
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FIG. 1: The classical autocorrelation function P{t) of Eq. is plotted for the perturbed sawtooth map of Eq. 114II . with 
A'l — 0.3, K2 = —0.2, and several values of v. P{t) is computed for each of the six observables B{q,p) — cos{2rmvq), sm{2mnq), 
where m = 1, 2, 3, and then averaged to obtain the behavior for a typical observable. The thin solid lines represent best fits 
to the long-time exponential behavior P{t) = foe"^"*, with A0.2 = 0.0083, A0.3 = 0.0275, A0.4 = 0.067, and A0.5 = 0.14. 



We may also consider a system with power- law decay of classical correlations, e.g., 

T{p) = -I- cos(27rp) 

f -i [92-^2/4] for|g|<7r/2 
I otherwise 

Here, bouncing-ball modes for |g| < tt/2 in the vicinity of p = dominate the long-time classical correlation decay, 
resulting in P{k) ^ bk~'^, with 7 = 1/3. The sum over time steps k in Eq. is effectively truncated on the scale 

of the Heisenberg time N by the smoothing function fe{k). Then J2k='i P{^)ff-{^) ~ N'^^'^ and thus 

-l(^)-7S- (19) 



IV. BROWNIAN MOTION MODEL FOR CHAOTIC ENERGY LEVELS 

The statistical properties of the energy levels £"„ can be obtained using a Brownian motion model P| . The discussion 
in this section follows |^ ^ |^ . In this model, the matrix elements of the Hamiltonian undergo a diffusive evolution as 
a function of a fictitious time variable r. We denote the infinitesimal change of H by 5H . The off-diagonal elements 
of 5H are assumed to satisfy 

&Hmn = , 6HmnSH*^,^, = 5t C°-^nS„n' 6mm' , (20) 

while the diagonal elements are taken to obey 

lff™:=0, SH^^SHnn ^^StC^^s^ (21) 

Within RMT, C^'^s = S^n and C°f„ = 1. Non-universal deviations from RMT are encoded in C^'^s = C^'!*; and 
^mn = ^m-n 13 ■ The Corresponding motion of the energy levels Eq , ... , E^^i may be analyzed as follows. Using 
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second order perturbation theory, one obtains 



5Hn 



(22) 



for the energy level shifts SEn- Thus 



^m—n 
En — Er, 



and 



■diag 
m — n ' 



(23) 



(24) 



It is convenient to express Eqs. H23() and H24|l in terms of the Fourier modes of Ai<^„ = En — nA. We enforce periodic 
boundary conditions Ai?jv — Ai?o and define 



N-l 



SO that 



N/2-1 



-Y ^En 



'2-nikn/N 



(25) 



n=0 



N/2-1 



2nikn/N , „ „ 

afc e ' + c.c. 



A 77 \ ^ 2'iTikn/N \ " 

i^En = ^ ake ' = } ^ 

k=-N/2 fe=0 

Since the A£'„ are real, = o,*-k- From Eq. (|24|l we have, using SAEn — SE. 



(26) 



6akSa* = SalSap ^ — St h Skp , 
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where h = J2n C^'^Se^'^^'^'"/^ , and 



SakSttp — 6al.Sa* — 



(27) 



(28) 



for k,p €z (0, N/2). In general, the expectation value of Sak is a complicated fmiction of all Op, which may be expanded 
as 



Sak — St 



Ai{k)ak + ^ A2{k, q)ak-qaq + ^ Az{k, q,p)aqapak-q-p + ■ 



(29) 



The coefficients Aj are obtained as follows (assuming that = 1, see however Ref. ^3)- Expanding the denom- 

inator En — Em in Eq. H23|) around the mean (n — m)A, 



SE„ 



^ iA 



1 - —{AEn+e ~ AEn 



2PA- 



■{AEn+t-AEnf 



(30) 



Again using SEn = SAEn and the identity e 



-2nikn/N 



X)?" ^k,jN, the first term in Eq. H3U|I implies 



Sak = ^Tat 



2-iTikl/N 



1 



(31) 



In the limit of large N one has to first order in |fc|/iV (extending the range of summation over ^ from — oo to oo) 

27r2|fc 



^ ^2^2 



7VA2 



(32) 
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and thus 

To lowest order in Eq. 129|l . the equilibrium distribution of the (corresponding to large fictitious time r) factorizes 
into a product of Gaussians, where the variance is given by 

2 NA^Ik 



while higher-order cumulants are zero. The approximation of Eq. H34|l is appropriate for k ^ N. In this case, the 
for different k are uncorrelated. In RMT, Ik = N^^, and 

2 A2 



for k<t: N. 

In order to determine the fluctuations of for fc ~ iV, higher-order terms in Eq. H29() must be taken into account. 
This may be done perturbatively, resulting in corrections to the variance and possibly non-zero higher-order cumulants. 
Moreover, for larger values of k may be correlated. 

In the above derivation, we have started with initial conditions ak = at fictitious time t — 0, i.e. we have assumed 
the initial unperturbed Hamiltonian has a "picket fence" spectrum. En = nA. A typical chaotic Hamiltonian, however, 
will not correspond naturally to a small perturbation of such a "picket fence" Hamiltonian. Therefore, in practice 
it is more useful to observe spectral equilibration by comparing coefhcients of a given Hamiltonian H with the 
coefficients a'j, of a perturbed Hamiltonian H + B. Full equilibration implies that the spectra of H and H + B become 
independent of one another, though drawn from the same random matrix ensemble. Then 



for /c < iV. 



V. CONDITIONS FOR ENERGY LEVEL EQUILIBRATION 

We now address more carefully the question originally raised in the seminal work by Wilkinson Q . Specifically, we 
wish to understand under what circumstances a class of classically small perturbations B of an initial Hamiltonian H 
is sufficient to generate random matrix statistics in the spectrum, at various energy scales. A related question is the 
size in parameter space of a random chaotic ensemble necessary to average away system-specific spectral properties 
and generate universal statistics 01 • 

Two conditions are necessary for equilibration to universal statistics to occur. First, if the perturbation is classically 
small, then equilibration to RMT may only occur on time scales longer than the decay time of classical correlations 
(otherwise, Ik for the corresponding modes a^). Thus, the modes ak may be equilibrated to RMT only for 

k ^ Tdocay/Tkick, independent of the perturbation, and correspondingly the spectrum may only display universal 
statistics on energy scales £ <C fi/T^cca.y Secondly, the perturbation must be sufficiently strong to equilibrate a given 
mode ak after the fictitious time r during which the Hamiltonian matrix elements undergo Brownian motion. In our 
units, this fictitious time r is simply equal to the variance 0: 

T = 4 , (37) 

and is given by Eq. (|10|l . On the other hand, the characteristic response time of the k-th mode is 

n _ 2h^ ^^g^ 

7rz^(i?)(fcTkick) (fc7kick)rH- ' 

as may easily be seen by comparing Eqs. 127|1 and (|34|l and then noting that the Heisenberg time Th is given by 
27r?i/A = 2TThiy{E) = NTi^ick- Equilibration will therefore happen for a given mode ak when t ^ Tfe, i.e. 



and transforming to the energy domain we find equilibration on scales 




(40) 



n 



The condition of Eq. H4UI) must be satisfied simultaneously with the first condition £ <^ ^/Tdocay, i-e. equilibration to 
universal statistics occurs on all energy scales £ satisfying 



Note that the two upper limits on the equilibration energy scale have opposite dependence on the classical correlation 

scale 2^dccay ■ 

Let us consider the behavior of the system as the size of the perturbation B is varied, for a given initial Hamiltonian 
H, at a given classical energy. We assume Tdocay < Th, so that the system is in the quantum chaotic regime, (i) For 
very small perturbations, {B"^) 7i^/(IdccayTi/), there is no equihbration on any time scale before the Heisenberg 
time Tfj and consequently no equilibration on any energy scale larger than a mean level spacing A. (ii) For larger 
perturbations, i.e. ^^(^decayT//) < (B^) < TiVT'locay aU modes k > fc^in - ?iV(Tkick7docay (B^)) equihbrate to 
their RMT values. In the energy domain, the corresponding condition is £ <^ ^equii (B')Tdocay/?i. (iii) Finally, for 
the largest (still classically small) perturbations fi^/T^ecay (-^^) ^ -^^i time scales beyond Tdecay and all energy 
scales below £oquii ^/Tdocay equilibrate to universal statistics. 

What happens if we instead fix the size of perturbation B, and consider a variety of classical systems, with different 
classical time scales Tdecay? (i) For the most strongly chaotic systems, Tdocay ^ h{B'^) , and equilibration is limited 
by the condition of Eq. (gU)). (ii) As the classical system becomes less chaotic, fcmin falls, more and more modes Uk 
come into equilibrium, and the energy scale £cquii increases, (iii) Maximum equilibration to RMT is attained when 

Tdecay ^{B^} ^^^i whcrc all modes Ofc for k > fcmin ^^kick(-^^) equihbrate, and £equii ~ (B^)^^^. (iv) Then, 
as the degree of chaoticity of the system continues to decrease, fcmin begins to increase, modes again move away from 
RMT equilibrium and the energy fequii begins to drop, (v) Eventually, we reach the border Tdecay ~ Th between 
quantum chaos and quantum regularity, where all equilibration to universal statistics is again absent. Thus, absence 
of equilibration to RMT at a given time or energy scale may be consistent with either very strongly chaotic or very 
weakly chaotic (or regular) systems, while maximum possible equilibration is attained for moderately chaotic systems 
between these two extremes. 

The discussion of the preceding paragraph is illustrated for the perturbed sawtooth map in Fig. [3 where we show 
the effect on the spectrum at various energy scales of a fixed-size perturbation, for four classical systems characterized 
by different decay times Tdecay The classical systems treated here are the same as in Fig. ^ while the family of 
perturbations used corresponds precisely to the family of classical functions used in that earlier figure. The squared 
change in each Fourier spectral coefficient ak has been normalized by the constant Cfj of Eq. so that full 

equilibration implies k\ak ~ a'y.\'^ /Cp = 1 for k <^ N. For k ^ N, this result is modified due to higher-order terms in 
Eq. H29I) . as discussed in Section Hvl and illustrated by the thin dotted curve in Fig.|2] From v = 0.70 to v = 0.20, the 
four systems clearly demonstrate the effect of gradually reducing the Lyapunov exponent and increasing the classical 
parameter Tdecay 

Looking first at the v = 0.70 curve in Fig. |21 we find equilibration of the spectrum only for k > 100, corresponding 
approximately to energy scales of three levels spacings or fewer. When the degree of chaos is reduced {v = 0.40) and 
Tdecay Correspondingly increases, the same perturbation strength is sufficient to equilibrate the spectrum up through 
much larger energy scales, namely k > kmin ~ 50. However, further reducing the degree of chaos by tuning v down 
below 0.4 and thereby increasing Tdecay serves to reduce the range of equilibrated energies for the same perturbation 
strength, as equilibration is now governed by the condition k ^ Tdecay/Tkick- Finally, for v = 0.20, Tdecay/Tkick — 120, 
and once again equilibration to universal behavior is not attained on any energy scale for our system size. 



In the previous section we considered spectral equilibration at various time scales fcTkick and energy scales E. We 
now focus specifically on the conditions for equilibration at the energy scale A, necessary for example to reproduce 
the Wigner-Dyson distribution of nearest neighbor level spacings. Assuming Tdecay < Th as before, we need to satisfy 
Eq. igni) for £ - A: 




(41) 



VI. EQUILIBRATION OF LEVEL SPACINGS 
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FIG. 2: The effect of a perturbation on the quantum spectrum as a function of Fourier mode k is plotted for four 
classical systems, distinguished by the parameter v of Eq. dJ. The values of Ki and K2 are as in Fig. Q and the 
Hilbert space dimension is A'^ = 512. As in the previous figure, the data are averaged over an ensemble of perturbations 
B{q,p) = X;l = , [a;™sin(27rmg) +j/™cos(27rmg)] with {B^) = 5 + 2^-) ^ O-GN-"" = 4.51 • 10'^ The constant Cf,, is 

given by Eq. 13611 . where (3 — 2 in the present case. The thin dotted curve indicates the limit of full equilibration on all energy 
scales, as obtained by choosing random and uncorrelated Hamiltonians H and H' — H + B . 



The size of the perturbation B must generically be at least of order to produce equiUbration of the level 

spacings. Less chaotic systems, however, as measured by a longer decay time, equilibrate more efficiently. For example, 
in a two-dimensional Hamiltonian system, B oc n^^'^ is needed for a fixed Tdocay, but B cx would be sufficient if 

Tdccay oc t]/^, and B cx h'^'"^'^ suffices if Tdccay oc Tjf. 

Equilibration at the scale of the mean level spacing may be measured by focusing on the Fourier coefficient a7v/2 
(notice that qn = ao)- Equilibration of this coefficient is studied in Fig.|31as a function of perturbation strength, for 
several classical Hamiltonians. Specifically, the mean squared fluctuation in the ajv/2 coefficient is plotted as a function 
of perturbation strength {B"^). For each classical system, we clearly observe the proportionality between |a7v/2 — a^^jP 
and (-B^) when |ajv/2— 1^/2!^ ^ N~^, and the eventual saturation at the system-independent equilibrium value at large 
perturbation size, just as predicted by the Brownian motion model. Furthermore, at a given perturbation strength, we 
see faster equilibration for systems with slower classical correlation falloff, as indicated by smaller classical parameter 
V. Semiclassically, we expect this increase in equilibration rate to be controlled by the integral dt P{t) « &Tdccay as 
in Eqs. ||SJ) and (|10|l : the classical predictions are indicated by dotted lines in Fig. 13 The slight discrepancy between 
the quantum data and the semiclassical prefactors may be explained by the finite system size, N = 256, since the 
classical expressions assume N ^ Tdecay/^kick, while in our case = 256 and Tdccay /Tkick reaches 120 when v = 0.2. 

Fig. ^ shows explicitly the dependence of equilibration rate on the classical system dynamics, by varying the 
classical parameter v while the perturbation strength is held fixed. Again, we see an order of magnitude change in 
the equilibration rate as v is varied, in agreement with the semiclassical prediction indicated by the dashed curve. 



VII. QUANTIZATION AMBIGUITY AND SPECTRAL EQUILIBRATION 

The original motivation behind Wilkinson's work 33 was to understand the relationship between spectral equi- 
libration and the quantization ambiguity, thus using the latter to explain the approach of generic chaotic systems 
to RMT behavior on short energy scales. Generally, quantization ambiguities come in two types. The first are 
gauge or boundary condition ambiguities that allow many different semiclassical theories to be associated with the 
same classical dynamics, as in the Aharonov-Bohm effect or in a change from Dirichlet to Neumann boundary con- 
ditions for a billiard system. These ambiguities correspond to 0(?ieff) terms in the quantum Hamiltonian, where 
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<b2> 



FIG. 3: Equilibration of the highest-frequency mode in the spectrum is studied for A'^ = 256 as a function of perturbation 
strength for several initial Hamiltonians. The classical systems as well as the ensemble of perturbations are the same as in 
Fig-El The thin dotted lines correspond to the classical prediction |ajv/2 ~ '^'n/2\ °^ (^^) /q dt P{t). 
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FIG. 4: Equilibration of the highest-frequency mode in the spectrum is studied, for TV = 1024 and perturbation strength 
(B^) = 6.4-10"^ as a function of the classical parameter v (solid curve). The classical systems and the ensemble of perturbations 
are the same as in Figs. Inland The classical prediction based on Eq. is indicated by the dashed curve. 

^off = ft/Styp ^ ?i/£'Tkick ~ f^/pL, Styp is the typical action of a short classical orbit, p is the typical momentum, and 
L is the system size. Secondly, there are operator ordering ambiguities that allow multiple quantum Hamiltonians 
to have the same semiclassical limit, including identical action phases at leading order. Ambiguities of this second 
class are 0(/i^ff) and result, for example, from canonically quantizing the same classical Hamiltonian in two coordi- 
nate systems related by a nonlinear canonical transformation 3]; they also appear naturally when different limiting 
procedures are used to define a quantum dynamics on a constrained surface |15|. 
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We first examine spectral equilibration due to the 0{h'^ff) quantization ambiguity, considered by Wilkinson in 
Ref. [3- Substituting B ~ ^cffE into Eq. (|41|l . we find equilibration on energy scales 

£ « fequii = £;min (n^^Ihi^Xs^) , (43) 

\ decay kick / 

or equivalently for Fourier modes 

k » k^,^ = max f K,'^) ■ (44) 

\ kick decay / 

For a fixed classical dynamics, the second expression in the parentheses in Eq. H43|) or Eq. 1)44(1 always dominates in 
the semiclassical limit TieS 0. For a d-dimensional Hamiltonian system, the total number of available fluctuating 
modes in the spectrum is 0{h^g'^), and all but the first 0(?i(J^Tkick/?dccay) of these are equilibrated. Thus, full 
equilibration is not possible for any spectral mode in the case d = 2, while for c? = 3 full equilibration is possible 
for the highest-fc modes, assuming slow decay of classical correlations, Tdccay/Tkick 3> 1. In the energy domain, this 
implies equilibration only on scales up to 0(Tdccay/2kick) mean level spacings A in the 3-dimensional case. The case 
d = 4 (e.g., two interacting particles in two dimensions) is the first for which equilibration generically extends to 
scales much larger than a mean level spacing, fequii 3> A. Finally, in the many-body limit d oo, an ever-increasing 
number of modes are equilibrated by the 0{1\.ff) ambiguity, however, the first 0(?i^Tkick/2dccay) modes are never 
equilibrated. Correspondingly the range of energies over which equilibration may occur always remains a factor of 
0(?ipffT'dccay/2kick) Smaller than the ballistic Thouless energy i^Thouicss ~ TicsE, independent of dimension. 

We now turn to the 0{fi~^) quantization ambiguity, associated with external gauge fields or boundary conditions. 
A similar analysis shows that equilibration now occurs on energy scales £ <C i^^cquii ~ ?ioff-E'2kick/7dccay and for Fourier 
modes k ^ fcrnin ~ Tdocay/Tkick- Thus, all modcs equilibrate in any dimension for a chaotic system, except for the 
first few that encode non-universal short-time dynamics for a classical system with slow correlation decay. 




FIG. 5: Equilibration of the highest-frequency mode in the spectrum of a perturbed sawtooth map is shown for 0{h^f[) 
perturbation strength (-B^) = (thick curves) and 0(?icff) perturbation strength (B^) — (0.16/A'^)^ (thin curves) as a 

function of A^. The thin dotted lines indicate the predicted A''"^ behavior for the 0(7icff ) case in the semiclassical limit A'^ oo. 

In Fig. 13 we again focus on equilibration of the highest-frequency Fourier mode ajv/2, corresponding to energy scales 
of order A, for our quantum map numerical model (which exhibits the scaling of a d = 2 autonomous Hamiltonian 
system). Consistently with the above discussion, an 0{1\.g) ambiguity is not sufficient to equilibrate even this highest- 
frequency Fourier mode in the semiclassical limit N oo (hcs 0); lower-frequency modes will be even further away 
from equilibration. On the other hand, an 0(?icff) (gauge-size) ambiguity provides full equilibration independent of 
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The above analysis and numerical simulation assume a constant Tdecay in the ?ioff — limit. However, a system 
remains in the quantum chaotic regime as long as Tdecay is shorter than the Heisenberg time, i.e. when Idecay/^kick < 
hl^"^. In general, the rate of equilibration to RMT behavior is increased by choosing a system with very slow classical 
relaxation, Tdecay ^ T^ick- For the d = 2 case, it is impossible to attain complete equilibration to RMT for any 
spectral mode with an 0{h^ff) perturbation, regardless of the choice of classical system. For d > 3, on the other 
hand, Eqs. H43|l and (|44|l clearly indicate that full equilibration for some modes can be achieved with an 0(?i^g) 
perturbation for systems with Tdccay/Tkick ^ 1 (where a generic chaotic system with Tdocay/Tkick ~ 1 would not 
display equilibrated behavior). Optimal equilibration is attained for Tdecay /Tkick ~ ^cS^ where all but the lowest 
0{h~g) modes are equilibrated in any dimension c? > 3. 
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FIG. 6: Equilibration of the highest- frequency mode in the spectrum is shown for O(ft^) perturbation strength (B^) = 
{0.66/ N^f, for a system with power-law classical correlation decay, defined by Eq. IjlSfl . solid curve, and for the perturbed 
sawtooth map with v = 0.7, dashed curve. The thin dotted lines indicate the predicted A^~^ behavior, for 7 = 1/3 and 'y — 1. 

Finally, in Fig. we compare equilibration in the hard chaotic perturbed sawtooth map with equilibration in a 
system governed by a slow power-law classical decay, defined by Eq. (|18|l . In both cases a perturbation of order 
t\.g ~ 7V~^ is used. Although no equilibration is possible for any system in d = 2, we see clearly that the power-law 
classical system comes closer to equilibration than the hard chaotic system for any large N , and the two systems are 
governed by different exponents in the — > 00 limit. 



VIII. FIDELITY AND SEMICLASSICAL ACCURACY 



The fidelity of quantum evolution in the presence of a small perturbation, also known as the Loschmidt echo, has 
received much attention recently, particularly in the context of classical-quantum correspondence and the very different 
behaviors that the quantum fidelity can exhibit for classically regular and chaotic systems For sufficiently large 
perturbations, the decay of quantum fidelity with time can be governed by the classical Lyapunov exponent 17, 'l^. 
Here we wish to mention an interesting connection between spectral equilibration and fidelity decay in quantum 
chaotic systems due to a small perturbation, as analyzed recently by Prosen and collaborators 9]. It was shown that 
for a static perturbation B with (B^) < h^/T^^^^^, the fidelity decays exponentially on a time scale varying inversely 

with Tdecay, namely Tfiddity ~ {B^)Tdccay Dividing through by the one-step time scale Tkick, we find that Tgdeiity 
corresponds to a dimensionless mode number 



(45) 
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which precisely agrees with the boundary between equihbrated and non-equiHbrated modes given by Eq. (|39(l . vahd 
as long as kadeuty > ^decay/Tkick- In other words, a mode k is equilibrated to universal behavior by a given class of 
perturbations if and only if two conditions hold simultaneously: the quantum fidelity has decayed to a value much 
less than unity by the corresponding time fcTkick cind this time scale fcTkick is larger than the classical relaxation time 

^decay ■ 

Another important connection is between spectral equilibration behavior discussed in the present work and the 
decay of semiclassical accuracy. Since different quantizations of the same semiclassical dynamics differ by 0(?ipg) in 
the Hamiltonian, the difference between a typical quantization and the semiclassical approximation must be at least 
of this order in the hes — > limit. Thus the error in the semiclassical approximation must be at least of the same size 
as the error caused by an 0(?i^ff ) perturbation in quantum mechanics, assuming of course that the physically correct 
quantization is "typical" . Clearly, hard quantum effects such as diffraction can increase the error in the semiclassical 
approximation, so the 0{t\.g) quantization ambiguity only provides a lower bound on the size of the semiclassical 
error, or equivalently an upper bound on the breakdown time of semiclassical validity. Previous work has shown that 
this bound on the semiclassical error is saturated for the case of a smooth Hamiltonian in the absence of caustics j^] , 
and furthermore it was demonstrated that the semiclassical error obeys different scaling laws with time and h for 
regular as opposed to chaotic classical systems, just as one would predict using the quantization ambiguity approach. 
The semiclassical accuracy problem may also be considered as an example of a generalized quantum fidelity problem, 
if non-Hermitian perturbations are considered (since scmiclassical evolution is in general non-unitary). 

The results of the present work, particularly those of Sec. IVIII imply an upper bound on the breakdown time 
scale of semiclassical accuracy: Tgomiciassicai ~ ^(Tff '^kick/^dccay and a lower bound on the breakdown energy scale of 
semiclassical accuracy: 

c pt3 ^decay p ^2 T^ccay 

C-scmiclassical ^ -t>"-cff "TTi -C'Thoulcss "-off "TTi : l^"7 

kick kick 

where -EThouioss is a ballistic Thouless energy. In particular, semiclassical accuracy will persist beyond the Heisenberg 
time Th ?i^7kick in any two-dimensional chaotic system (assuming diffraction and caustics are properly accounted 
for), allowing individual energy levels and wave functions to be semiclassically resolved. In three-dimensional systems, 
Th ^ ?icff ^kick, and Tscmidassicai ~ T/flkick/Tdocay ■ Here, the semiclassical breakdown time may be comparable to the 
Heisenberg time for the most chaotic systems (Tdecay/Tkick ~ 1), but any slowdown in the classical relaxation rate will 
cause a faster breakdown in the semiclassical approximation. Finally, given four or more classical degrees of freedom, 
semiclassical accuracy inevitably breaks down well before the Heisenberg time for any generic dynamics, even the 
most chaotic, and semiclassical reproduction of individual spectral levels is never possible. It may be of interest to 
investigate the manner in which improved semiclassical approximations beyond leading order in h |19| may produce 
different scaling of the semiclassical error and possibly permit quantum spectra to be semiclassically resolved in four 
dimensions and higher. 



IX. SUMMARY 



A careful examination of the Brownian motion model for spectral equilibration to universal statistics shows that 
equilibration is strongly dependent on the classical relaxation rate as well as on the dimensionality of the system. Two 
competing time scales, proportional and inversely proportional to the classical relaxation time Tdocay, jointly govern 
the equilibration process. Balancing of these two time scales implies that for a given perturbation B, equilibration of 
the maximum number of spectral modes is achieved when Tdccay h{B^)~^/'^ . For small perturbations, B < h/Tdccay, 
a relation exists between spectral equilibration as a function of mode number k and the decay of the quantum fidelity 
as a function of time. 

Focusing on the effect of O^h^g) perturbations, associated with the ambiguity of quantization, we find that no 
equilibration to universal statistics is ever possible for dynamics in two dimensions. In three dimensions, equilibration 
to universal statistics occurs for some modes, but only if the classical relaxation time Tdccay is sufficiently long, while 
in four dimensions and higher, equilibration of at least some modes occurs for any chaotic system. However, optimal 
equilibration only occurs for very long classical relaxation times, Tdccay oc h~g. These predictions of the Brownian 
motion model are also entirely consistent with results for semiclassical accuracy in smooth chaotic systems. 
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